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If a galaxy resides in a cluster, then its passage through the pervasive intracluster 
medium will produce a detectable signature in the X-ray emission from the cluster. 
Such features have now been detected in a number of systems. The simplest kinematic 
information that can be extracted from this signature is the galaxy's direction of 
motion on the plane of the sky. This paper explores the constraints on cluster dynamics 
that could be derived from such information. In particular, we show that it is possible 
to define a projected anisotropy parameter, B(R), which is directly analogous to the 
usual orbital anisotropy parameter. We describe an estimator for this quantity, B(R), 
which can be derived in a robust and straightforward manner. We present a simple 
dynamical model for a cluster consisting of a Michie distribution function of galaxies 
orbiting in a truncated singular isothermal sphere potential. Using this model, we 
demonstrate the ambiguity between the distribution of mass and the distribution 
of galaxy orbits when interpreting the traditional measures of cluster kinematics (the 
projected density of galaxies and their line-of-sight velocity dispersion). As an example, 
we show how two very different dynamical models can fit the kinematic properties of 
the Coma cluster. We demonstrate that the measurement of B using a relatively 
small sample of wake directions (N ws ^ c w 50) would provide an effective mechanism 
for lifting this degeneracy. Thus, by combining X-ray measurements of wake directions 
with number counts and line-of-sight velocities derived from optical data, it will prove 
possible to measure both the orbit distribution and the form of the gravitational 
potential in clusters of galaxies. The requisite X-ray observations lie within reach of 
the soon-to-be-launched AXAF satellite. 
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1 INTRODUCTION 

As the seminal work on the Coma cluster by Kent & Gunn 
(1982) first demonstrated, it is possible to study the dy- 
namics of a cluster of galaxies by combining photometric 
observations of the galaxies' positions on the plane of the 
sky with spectroscopic measurements that reveal the line- 
of-sight velocities of cluster members. If it is assumed that 
the cluster is spherically symmetric, these observed prop- 
erties can usefully be quantified by the projected number 
density of galaxies as a function of projected radius, N(R), 
and the RMS dispersion in the galaxies' line-of-sight veloc- 
ities as a function of projected radius, a\ OB (R). Dynamical 
models with different distributions of galaxy orbits can then 
be tested against the observations by comparing N(R) and 
o~i os (R) to the forms for these functions that the models pre- 
dict. 



In principle, such a study allows us to investigate the 
manner in which a cluster has formed and evolved. Galax- 
ies have relatively long dynamical "memories," and so their 
current orbits must to some degree reflect the processes by 
which the cluster formed. For example, if the galaxies are 
found to follow preferentially radial orbits, then their mo- 
tions presumably reflect the radial infall from which the clus- 
ter originally formed. Further, the kinematics of the galaxies 
are dictated by the gravitational potential of the cluster as a 
function of radius, $(r), and so dynamical modelling allows 
us to probe the distribution of mass in the cluster. Since the 
mass of a galaxy cluster is dominated by dark matter, no 
study of the evolution of such a system is complete with- 
out addressing the distribution of this massive component 
as inferred from the gravitational potential. 

Unfortunately, there is a fundamental ambiguity in such 
studies. As Binney & Mamon (1982) demonstrated, for any 
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given <6(r) it is possible to reproduce a wide range of func- 
tions N(R) and ai oa {R) simply by varying the balance be- 
tween radial and tangential motions in the orbits that the 
galaxies follow. They parameterized this balance by the 
anisotropy parameter, 

/3(r) = l-^J, (1) 

where at and ay are the RMS velocities in the tangential 
and radial directions for galaxies at radius r in the cluster.^ 
Thus, it is not possible to use just N(R) and ai os (R) to solve 
for both $(r) and /9(r). 

One approach to lifting this degeneracy is to find some 
independent technique for measuring <&(r), and then use 
N(R) and ai OB (R) to solve for the orbital anisotropy, /3(r). 
For example, the form of the gravitational potential can be 
inferred from the properties of the X-ray emitting intraclus- 
ter medium (ICM). By measuring the temperature of this 
hot gas and its spatial distribution using an X-ray telescope, 
it is possible to solve uniquely for the form of the gravita- 
tional potential that contains the gas - see, for example, 
Eyles et al. (1991). Alternatively, one can use the distortion 
of optical images of background galaxies due to the gravi- 
tational lensing effects of the cluster mass to measure <I>(r) 
directly (Tyson, Wenk & Valdes 1990). 

In practice, neither of these techniques is entirely satis- 
factory. The measurement of $(r) using X-ray data requires 
one to measure the temperature profile of the ICM as well 
as its density profile, and the poor spatial and spectral reso- 
lution and limited sensitivity of existing X-ray observations 
make this a difficult task. The inversion from gravitational 
lens distortions to 3>(r) can only be made if the distorted 
shapes of many background galaxies are measured, and is 
to some extent dependent on the assumed redshift distri- 
bution of the background galaxies. Further, neither of these 
techniques is well-suited to measuring &(R) in poorer clus- 
ters, where the X-ray emission will be weak, and the amount 
of gravitational-lens distortion will be small. 

Perhaps this problem can be more profitably attacked 
from the other direction: if we can find some additional con- 
straint on the distribution of orbits in a cluster, then we 
can use this information to lift the ambiguity and solve 
for both the galaxy kinematics and the cluster potential. 
One attempt to adopt this approach was made by Merritt 
(1987). In a further study of the Coma cluster's kinematics, 
he showed that the distribution of member galaxies' line- 
of-sight velocities will have a shape that may depart dra- 
matically from a Gaussian depending on the distribution of 
galaxy orbits. Specifically, if the orbits are preferentially ra- 
dial then the the velocity distribution will have longer tails 
than a Gaussian, while if the orbits are more circular it will 
have a "boxier" appearance. This information in the shape 
of the line-of-sight velocity distribution has been success- 
fully exploited in dynamical studies of the stars in individ- 
ual galaxies (e.g. van der Marel & Franx 1993, Kuijken & 
Merrifield 1993) , but it has proved less useful in the case of 
clusters. The main difficulty in using the shape of a cluster's 

* Note that there are two equal tangential components in a spher- 
ical system (in the 9 and </> directions of a spherical polar coor- 
dinate system), and so cr^ is sometimes defined to be a factor of 
two larger. 



velocity distribution is that clusters are not entirely relaxed 
systems. A relatively small amount of substructure in the 
distribution of galaxy velocities within a cluster, perhaps 
due to an infalling group of galaxies that has yet to dis- 
rupt completely, can totally dominate the departures from 
a Gaussian in the line-of-sight velocity distribution (Zablud- 
off, Franx & Geller 1993). 

In this paper, we explore an alternative candidate for 
the requisite kinematic constraint: the X-ray signature of 
the interactions between cluster galaxies and the surround- 
ing ICM. As a galaxy moves through the ICM, the combina- 
tion of gravitational focusing of the ICM and ram pressure 
stripping of the galaxy's own interstellar medium (ISM) will 
produce an enhancement in the density of gas behind the 
galaxy - see Balsara, Livio & O'Dea (1994) for a review of 
these processes and the results of some simulations. Thus, 
we expect a galaxy's direction of motion on the plane of the 
sky to be marked out in X-ray observations by a "wake" 
of excess emission trailing behind it. This phenomenon has 
now been observed in several systems. Sakelliou, Merrifield 
& McHardy (1997) found a wake of X-ray emission emanat- 
ing from the radio galaxy 4C34.16. In this case, the bent 
morphology of the radio jets provides a further measure of 
the galaxy's direction of travel, and the wake is found to lie 
downstream from the moving galaxy. A further example is 
provided by the dumb-bell galaxy NGC 4782/3 (Colina & 
Borne 1995). In this interacting pair of elliptical galaxies, 
the X-ray emission from the individual components is offset 
from the optical light produced by the galaxies in the sense 
that one would expect if it were trailing behind the orbit- 
ing galaxies. The X-ray emission also reveals more extended 
wake features lagging behind the orbiting components. Fi- 
nally, Jones et al. (1997) have found that the X-ray emission 
from the elliptical galaxy NGC 1404 is strongly distorted, 
forming a wake pointing away from the centre of the For- 
nax cluster in which it resides. They suggest that this X-ray 
morphology arises from the ram pressure stripping of the 
galaxy's ISM, and hence conclude that NGC 1404 is on a 
plunging radial orbit. 

With the advent of more sensitive X-ray telescopes such 
as XMM and high resolution imaging instruments such as 
AXAF, the detection of these wake phenomena can only 
become more commonplace, and we might reasonably expect 
to be able to use wakes as a probe of the kinematics of many 
galaxies within a single cluster. For example, Balsara et al. 
(1994) have calculated that at the distance of the Coma 
cluster one might expect to see many wake structures with 
extents of ~ 2 arcseconds - a scale readily resolvable with 
AXAF. 

In principle, a great deal of information can be gleaned 
from the pattern of shocks and other structure present in 
X-ray wakes. However, in addition to the velocity of the 
galaxy, this structure depends in a complex manner on the 
properties of the ICM, how much of the galaxy's ISM has 
already been stripped away, how rapidly the ISM is being 
replenished by stellar mass loss within the galaxy, etc. We 
therefore content ourselves with using the X-ray morphol- 
ogy to measure the direction in which each galaxy is moving 
on the plane of the sky. In the remainder of this paper, 
we investigate whether such information is sufficient to lift 
the degeneracy between orbital anisotropy and gravitational 
potential. In Section 2, we present a practical measure of 
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the orbital structure that can be determined directly from 
the observed wake directions in a cluster. Section 3 presents 
a simple dynamical model, which illustrates how the pro- 
posed measure of orbital structure can discriminate between 
otherwise-indistinguishable models. 



2 ESTIMATING THE PROJECTED 
ANISOTROPY 

The direction in which the wakes in a cluster trail on the 
sky clearly tells us something about the distribution of or- 
bits. Suppose, for example, that the wakes imply that all 
the galaxies in the cluster are moving either directly toward 
or away from the cluster centre in the plane of the sky. The 
only spherically-symmetric system that could produce such 
an arrangement of projected motions is one in which the 
galaxies are all on intrinsically radial orbits. Thus, in this 
case at least, the distribution of wake directions is sufficient 
to define the distribution of orbits uniquely. Note, however, 
that the relation between the wake directions and orbital 
anisotropy is not completely trivial: if all the galaxies follow 
circular orbits, the wakes will not all be tangentially-oriented 
on the plane of the sky, and some will still appear pointed 
at the cluster centre. 

In order to explore the relationship between wake direc- 
tions and orbital distribution more fully, it is useful to define 
a "projected anisotropy." By analogy with equation ([l]), an 
obvious definition for this quantity is 

B(R) = 1-4- (2) 

where ctt(-R) and or(R) are the RMS velocities of galaxies 
on the plane of the sky in the projected tangential and radial 
directions, respectively. 

The projected anisotropy cannot be measured directly, 
but we can estimate it from the observed distribution of pro- 
jected orbit directions. If we define O to be the angle between 
the direction of a galaxy's motion on the plane of the sky 
(as inferred from its wake) and the projected radius vector 
of the galaxy, then the quantity that we can observe directly 
is the distribution of projected orbit directions as a function 
of projected radius, P(Q,R). This quantity can be related 
directly to B(R). For example, let us assume that at some 
projected radius in the cluster, the plane-of-sky velocity dis- 
tribution (POSVD) can be approximated by a function of 
the general form 

F(vr, vt) — F(v R /a 2 + v T /b 2 ), (3) 

where a and b are constants. By transforming from the coor- 
dinates {vr, vt} to polar coordinates on the plane of the sky, 
{i> pos , O}, it is straightforward to show that the distribution 
of observed values of O will be 



P ^ = 2,(1 -iW 8) - (4) 

Thus, we can estimate B by fitting the functional form of 
equation (Q) to the observed distribution for P(O) at any 
given projected radius. 

One robust method for implementing this fitting pro- 
cedure is to measure the fraction of galaxies for which 
7r/4 < 9 < 37r/4 or 57r/4 < 6 < 77r/4. This quantity, 



/t, measures the fraction of galaxies whose orbits are closer 
to tangential than radial on the plane of the sky. By inte- 
grating equation (^), it is straightforward to show that the 
value of B that gives rise to a particular value of fx is 

B = l-tan 2 (|/ T ). (5) 

Thus, by measuring fx, we can use equation (^) to estimate 
the projected anisotropy. This estimator for B has a number 
of desirable properties: it can be readily calculated from ob- 
served wake directions; it is robust in the sense that a small 
fraction of badly-determined directions will not affect it un- 
duly; and it has error properties which can be determined 
straightforwardly from the binomial distribution of fr- 

Strictly speaking, B is only an unbiased estimator for 
B if the POSVD takes the elliptical form of equation (^J). 
However, a rather wide range of plausible models can be ap- 
proximated in this way (see below). Even in cases where this 
approximation is not valid, B still provides a quantity that 
can be readily determined not only for observational data, 
but also for dynamical models containing different orbit dis- 
tributions. Thus, by comparing the observed profile B(R) to 
that predicted by a dynamical model, it is possible to test 
the model against observation. 



3 CASE STUDY: THE MICHIE/TRUNCATED 
SINGULAR ISOTHERMAL SPHERE MODEL 

To see how useful B(R) might prove as a diagnostic of orbital 
structure, we now consider a simple family of dynamical 
models. Such models are defined by the distribution function 
for the cluster, /, which specifies the density of galaxies as a 
function of both velocity and position (i.e. the phase density) 
- see Binney & Tremaine (1987) for a thorough discussion of 
distribution functions. By Jeans theorem, for a spherical sys- 
tem in equilibrium, the distribution function depends only 
on the energy, E, and angular momentum, J, of the point 
in phase space under consideration. A relatively simple dis- 
tribution function that might approximate the properties of 
a real cluster of galaxies is the Michie model, 

m j) = ( const.[e-^ - l]e- j2 ^ 2 o^ E<0 (g) 
[_ otherwise, 

where oq and r a are constants (Michie 1963, Binney & 
Tremaine 1987). This function produces a model where the 
distribution of orbits is isotropic at small radii, but becomes 
radial at large radii. Thus, it mimics what we might expect 
in a cluster of galaxies: relaxation effects will have random- 
ized the orbits at small radii in the cluster, but galaxies that 
travel to large radii will not have relaxed, and so their orbits 
will still reflect the radial infall by which the cluster formed. 
The quantity r a defines the characteristic radius at which 
the transition from isotropic to radial orbits occurs. 

The other dynamical quantity that we must specify for 
this model cluster is the gravitational potential, &(r), by 
which the galaxies are confined. Since clusters of galaxies 
are dominated by dark matter, and since we are interested 
in varying the orbital anisotropy and gravitational potential 
independently, we do not make the customary assumption 
that mass follows light. Instead, we impose an external grav- 
itational potential, and vary its parameters to assess the im- 
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Figure 1. The observable kinematics for the dynamical mod- 
els discussed in the text (with parameters ro and r a as anno- 
tated). Plotted as functions of projected radii, the upper panel 
shows projected density, and the middle panel shows the line- 
of-sight velocity dispersion. The lower panel shows the projected 
anisotropy: the lines show the estimator calculated using equa- 
tion (^|), while the crosses mark the corresponding true projected 
anisotropy [equation (^)] . The central projected densities and line- 
of-sight velocity dispersions have been scaled to a common value. 



pact of changing the distribution of mass in the cluster. The 
simplest potential that enables us to carry out such an inves- 
tigation is the truncated singular isothermal sphere (TSIS), 
which is produced by a mass density distribution 



p(r) 



a 2 /(2irGr 2 ) 




r < r 
otherwise, 



(7) 



where the constant ao sets a characteristic velocity scale for 
this potentialjj] The TSIS density distribution produces a 
gravitational potential 



$(r) = 



2 ^o [In (£) + - 1] r<r 
-2alr (i - ^-5— ) otherwise. 



(«) 



Here, we have set the arbitrary constant offset in the poten- 
tial such that $(r > r max ) > 0. Thus, from equation (^), all 
the members of the cluster in this model are confined to the 
region r < r max . 

Given / and we can calculate the density and velocity 
distribution of galaxies at any point in the cluster. It is thus 
straightforward to calculate observable quantities - such as 
N(R), £Ti os (J?), and the projected anisotropy estimator B(R) 
- by integrating the properties of the cluster along lines 
of sight through the system at different projected radii R. 
Figure |l| shows the results for a set of such calculations for 



t There is no fundamental reason why this characteristic velocity 
should be the same as the one in equation (Q), but we adopt this 
value as the simplest choice. 



the Michie distribution function in the TSIS potential. The 
spatial scale has been set by defining r max = 1, and various 
combinations of the mass cut-off radius, ro, and the orbital 
anisotropy radius, r a , have been adopted. 

For the case {ro = 1, r a — > oo}, the mass cut-off has no 
impact (since it lies outside the cluster), and so the poten- 
tial is that of a singular isothermal sphere. At small radii, 
the Michie distribution function is indistinguishable from 
an isothermal distribution function, and so the cluster itself 
mimics a singular isothermal sphere, with N oc R~ and ai os 
approximately constant. At larger radii, the more energetic 
galaxies, which would have been included in an isothermal 
distribution, are excluded from the Michie model, and so N 
drops faster than and <ri oa starts to fall. Both these 

quantities reach zero at R = r max , where the cluster ends. 

As ro is decreased, the truncation of the mass distri- 
bution affects both TV and <ti os . The models with ro = 0.3 
and ro = 0.15 have insufficient mass to contain the highest- 
velocity cluster members, and so the density profile and 
the velocity dispersion both drop more rapidly with radius. 
Thus, both the dispersion profile and density profile decline 
rapidly when the edge of the mass distribution is reached. 

However, as discussed in the introduction, the observed 
kinematics are affected by orbital anisotropy as well as the 
shape of the potential. As Fig. [j] illustrates, when r a is re- 
duced to a finite value, the absence of galaxies on tangen- 
tial orbits at large radii means that the projected density 
drops more rapidly than for the isotropic case. Further, since 
galaxies on radial orbits travel mostly transverse to the line- 
of-sight when seen at large projected radii, the observed line- 
of-sight velocity dispersion is also depressed. 

In fact, as Fig. |l| shows, it is possible to produce very 
similar profiles for N(R) and ai os (R) using very different 
combinations of ro and r a . In the case of the {ro = 1.0, r a = 
0.1} and the {ro = 0.15, r a — > oo} models, for example, the 
values of N differ by less than ~ 10% inside R = 0.8r max . 
Within this radius, the values of <ti os returned by the two 
models differ by less than 5%. Models with less extreme 
choices of these parameters, which include both a cut-off in 
the density profile and some degree of orbital anisotropy, 
differ even less. 

The futility of trying to distinguish between such mod- 
els is illustrated in Fig. |^. This figure shows the kine- 
matic properties of the Coma cluster as derived by Kent 
& Gunn (1982) using line-of-sight velocities of more than 
200 galaxies, with two of the models from Fig. |l] super- 
imposed. It should be stressed that these models have not 
been formally fitted to the data, but rather they are simply 
the arbitrarily-chosen models calculated above, with their 
number-, velocity-, and length-scales chosen to match the 
observations. Indeed, it would appear that the models fail 
to reproduce the observed flattening of the density profile at 
small radii. However, Beers & Tonry (1986) have shown that 
such a flat core is an artifact that arises from errors in the 
choice of cluster centre. With the centres measured correctly, 
they discovered that the projected density profiles of most 
clusters contain no core at small radii, but continue to rise 
in a power-law cusp, TV oc -R -1 , just as the models do. What 
is clear from Fig. |2| is that there is no way that the very dif- 
ferent dynamical models can be discriminated between on 
the basis of these data. With very much larger data sets, 
the error bars could be beaten down to a point where the 
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log(i?/arcmin) 

Figure 2. The projected density and line-of-sight velocity dis- 
persion as a function of projected radius for Coma cluster galaxies 
(Kent & Gunn 1982) as compared to two of the models discussed 
in the text. 



differences might become significant. However, it is already 
apparent from the lower panel of Fig. ^ that there are likely 
to be systematic errors in the observed properties: modest 
amounts of substructure in the cluster will mean that the 
simple spherically-symmetric model is at some level invalid, 
and contamination by unrelated galaxies along the line of 
sight may bias the derived properties. The presence of these 
fundamental systematic uncertainties means that the collec- 
tion of the vast data set necessary to beat down the random 
errors would be a nugatory exercise. 

Figure ^ shows how the measurements of wake direc- 
tions, and hence B(R), might be used to lift this effec- 
tive degeneracy between the different dynamical models. 
For the models in which r a — » 00, the distribution of or- 
bits is isotropic at every point in the cluster. Thus, in these 
cases there are as many galaxies moving on radial orbits as 
tangential ones on the plane of the sky, the distribution of 
wake directions is isotropic, and hence B = 0. However, for 
the models with finite values of r a , the radial orientation 
of the orbits is directly reflected by the positive values of 
the projected anisotropy parameter. Although the POSVD 
in this model is not of the form specified by equation (|^), 
it is clear that B(R) is nonetheless a good approximation 
to B(R) for this moderately-realistic distribution function; 
the two quantities never differ by more than ~ 10%. Thus, 
the observable B(R) does provide a useful measure of the 
underlying projected anisotropy. 

The finite values of B predicted by the anisotropic mod- 
els would be readily detectable with relatively small sam- 
ples of wake directions. If, for example, we assume that 
the cluster is isotropic, then fr = 0.5, and we can apply 
binomial statistics to show that a sample of N^ake mea- 
sured wake directions will produce an estimate for this quan- 
tity of f T = 0.5 ± 0.5/VN wa „ kc . The value of B « 0.7 in 



the anisotropic models shown in Figure |l| corresponds to 
fr ~ 0.3. Thus, if the cluster were in reality better mod- 
eled by one of these anisotropic systems, then a sample of 
only -/Vwake w 50 wake directions drawn from projected radii 
where B w 0.7 would be sufficient to rule out the isotropic 
fr = 0.5 model at the 3<r level. 



4 DISCUSSION 

The passage of a galaxy through its host cluster must have 
some impact on the distribution of the gas in its vicinity, 
and we are now finding the first signs of these effects in 
X-ray observations. The purpose of this paper has been to 
address the question of whether the simplest kinematic in- 
formation that can be extracted from observations of this 
phenomenon - the distribution of directions of galaxy mo- 
tion on the plane of the sky - can be useful in dynami- 
cal studies of clusters. By means of the simple Michie/TSIS 
model, we have shown that such observations can readily 
discriminate between models that have different underlying 
orbit distributions and gravitational potentials, but that are 
indistinguishable using traditional kinematic analyses. 

In practice, the proposed measure of orbital anisotropy 
on the plane of the sky, B(R), may not make optimal use of 
the observations, since it requires that we bin the data radi- 
ally. It also discards any information that might be gleaned 
from coupling between directions on the plane of the sky 
and line-of-sight velocities: it is, for example, possible that 
galaxies with the largest velocities might lie preferentially on 
radial orbits, and so there might be a positive correlation be- 
tween line-of-sight velocity and how close a galaxy's orbit is 
to radial on the plane of the sky. It would be straightforward 
to exploit this information by doing a full maximum likeli- 
hood fit to the observed wake directions and line-of-sight 
velocities of a sample of galaxies in a cluster. The advantage 
of the approach adopted here is that B can be calculated in 
a manner that is robust against the effects of a few bad data 
points, while such erroneous data could seriously compro- 
mise a maximum likelihood analysis. Further, the close tie 
between B(R) and the intrinsic anisotropy parameter, /3(r), 
means that this simple-to-measure quantity can readily be 
interpreted in terms of the orbital properties of the cluster. 

It is interesting to note that the signature of anisotropy 
in a cluster is likely to be found at large radii (see Fig. hi). 
This radial dependence is fortunate, since it will be easiest 
to detect a wake near the outskirts of the cluster, where 
we are not viewing the interactions between the galaxy and 
the ICM through a large column of unrelated X-ray emit- 
ting gas. Similarly, this technique is particularly appropri- 
ate to poor systems, where the column of ICM along the 
line of sight is relatively small. Further, in poorer systems 
the energy density stored in the ICM is lower, and so the 
disturbance caused by the passage of a galaxy will be larger 
and hence more readily detectable. It is also noteworthy that 
contamination by non-cluster galaxies along the line of sight 
is one of the major problems that dogs studies of low den- 
sity regions such as poor clusters and the outskirts of rich 
clusters. Since such unrelated galaxies will not interact with 
the ICM, it will be possible to screen them out from the pro- 
posed studies based on their lack of wakes. The applicability 
of the technique to large radii in rich clusters and to poor 
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clusters is particularly exciting, since, as was discussed in 
the introduction, it is in exactly these areas that the exist- 
ing methods for mapping the gravitational potential break 
down. 

One possible difficulty in using wake directions to study 
cluster dynamics is that the detectability of a wake may vary 
systematically with the dynamics of the galaxy producing it. 
For example, galaxies traveling rapidly through the cluster 
will produce longer and more dramatic wake phenomena, 
which will be easier to observe. Similarly, galaxies that have 
only recently joined the cluster are more likely to contain a 
significant ISM, which will be readily detectable as a wake 
of stripped material. This possible bias can be minimized 
by obtaining as deep an X-ray image as possible, so that 
even the fainter wakes are detected. Fortunately, an imag- 
ing X-ray telescope will record many wakes simultaneously, 
and so the mapping process can be carried out efficiently 
in few exposures. Hopefully, with its high spatial resolution 
and sensitivity, AXAF will open this new window on the 
dynamics of clusters of galaxies. 



ACKNOWLEDGEMENTS 

The author gratefully acknowledges the support of a PPARC 
Advanced Fellowship (B/94/AF/I840). 



REFERENCES 

Balsara, D., Livio, M. & O'Dca, CP., 1994, ApJ, 487, 83 

Beers, T.C. & Tonry, J.L., 1986, ApJ, 300, 557 

Binney, J. & Mamon, G.A., 1982, MNRAS, 200, 361 

Binncy, J. & Trcmainc, S., 1987, Galactic Dynamics. Princeton 

University Press, New Jersey 
Colina, L. & Borne, K.D., 1995, ApJ, 454, L101 
Eyles, C.J., ct al., 1991, ApJ, 376, 23 
Jones, C, et al., 1997, ApJ, 482, 143 
Kent, S.M. & Gunn, J.E., 1982, AJ, 87, 945 
Kuijken, K. & Merrifield, M.R., 1993, MNRAS, 264, 712 
Merritt, D., 1987, ApJ, 313, 121 
Michie, R.W., 1963, MNRAS, 125, 127 

Sakelliou I., Merrifield M. R., McHardy I. M., 1996, MNRAS, 283, 
673 

Tyson, J.A., Wenk, R.A. & Valdcz, F., 1990, ApJ, 349, LI 
van der Marel, R.P. & Franx, M., 1993, ApJ, 497, 525 
Zabludoff, A.I., Franx, M. & Geller, M.J., 1993, ApJ, 419, 47 



